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We discuss a method for scaling a neutral atom Rydberg gate quantum processor to a large number 
of qubits. Limits are derived showing that the number of qubits that can be directly connected by 
entangling gates with errors at the 10 -3 level using long range Rydberg interactions between sites 
in an optical lattice, without mechanical motion or swap chains, is about 500 in two dimensions and 
7500 in three dimensions. A scaling factor of 60 at a smaller number of sites can be obtained using 
collective register encoding in the hyperfine ground states of the rare earth atom Holmium. We 
present a detailed analysis of operation of the 60 qubit register in Holmium. Combining a lattice 
of multi-qubit ensembles with collective encoding results in a feasible design for a 1000 qubit fully 
connected quantum processor. 
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I. INTRODUCTION 

Quantum computing has the potential for perform- 
ing numerical calculations such as factoring and un- 
structured search faster than is possible on classical 
computers Despite the slower speed of basic gate 
operations on quantum compared to classical computers 
an algorithmic speedup is predicted for sufficiently large 
problems. Impressive yet rudimentary demonstrations 
of small quantum algorithms using less than 10 qubits 
have been achieved using several different physical em- 
bodiments of quantum bits and gates. These include 
controlled entanglement of up to 8 trapped ionsJH, 
and factoring with molecular spins manipulated by nu- 
clear magnetic resonance techniques [J]. In order to real- 
ize the full potential of quantum algorithms for solving 
difficult computational problems it will be necessary to 
develop approaches which allow a large number of qubits 
to be interconnected. Several conceptual designs for scal- 
able quantum computing architectures have appeared in 
recent years based on ion traps [5|, and there is intense 
effort directed at scalability of other approaches includ- 
ing quantum dotsmfl, superconductors [7[, linear opti cs [8j , 
rare earth crystals [9|, and small quantum repeaters [10(. 

In this paper we examine the potential scalability of a 
computer based on neutral atom qubits trapped in op- 
tical lattices with long range two-qubit gates mediated 
by dipolar interactions of Rydberg states[ll|. We will 
carefully estimate the maximum number of qubits that 
can be directly interconnected. Although it is in princi- 
ple possible to move information arbitrary distances in 
a quantum computer by swap operations between neigh- 
boring qubits, there is strong evidence for an increase in 
the error threshold compatible with scalable computa- 
tions when swap operations are required [12j. Conversely 
it is possible to put rigorous bounds on the required error 
threshold in globally connected models An alterna- 



tive approach to coupling qubits separated by distances 
that exceed the range of a direct interaction is based 
on conversion of the quantum information from station- 
ary to flying and then back to stationary qubits. This 
potentially eliminates the need for a long chain of swap 
operations but introduces the nontrivial requirement of a 
high efficiency stationary - flying qubit interface for most 
of the qubits in the quantum processor. 

While we do not dispute the potential of the above ap- 
proaches to long range interconnections it is interesting 
to explore how large a system one might build that allows 
each pair of qubits to be directly entangled without inter- 
mediate steps. We present here an analysis that takes ad- 
vantage of the long range of Rydberg atom dipole-dipole 
interactions as well as the scaling factor provided by col- 
lective addressing [LH Il5j. We divide our discussion into 
several parts. In Sec. |TT] we make estimates of the num- 
ber of qubits that can be directly connected in an optical 
lattice using Rydberg interactions. In Sec. IIIII we show 
that collective encoding in the rare earth atom Ho can 
potentially provide a scaling factor as large as 60 per site. 
In Sec. IIVI we discuss a protocol for filling an array of 
localized ensembles in a planar lattice, and estimate the 
total size of the globally connected processor. We con- 
clude in Sec. [V] with a summary and outlook for the 
future. 



II. RYDBERG GATES IN OPTICAL LATTICES 

The idea of using the dipole-dipole interaction of Ry- 
dberg atoms for neutral atom quantum gates was intro- 
duced by Jaksch and coworkers in [Tl| and subsequently 
extended to a mesoscopic encoding of the qubit by Lukin 
and coworkers p^j. A number of subsequent papers have 
analyzed in more deta il sp ecific schemes for implement- 
ing Rydberg gates [13, EH Ell with the conclusion that 
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FIG. 1: (color online) Planar optical lattice denning an array 
of sites with spacing D, each one of which may be a sin- 
gle atom or a small ensemble containing K atoms. Rydberg 
coupling enables gates within each ensemble and betweeen 
ensembles. 



they present a very promising approach to quantum logic. 
At this time a neutral atom Rydberg gate has not been 
demonstrated experimentally, although many of the un- 
derlying requirements have been separately achieved in- 
cluding loading and manipulation of single atoms in op- 
tical traps [2fJ,|2jJ,|22|, signatures of Rydberg interactions 
and dipole blockade at the many atom level !23i], and 
recently coherent excitation of Rydberg states together 
with observation of interaction effects at the level of two 
atoms [IH. 

We envision an array of atoms in optical traps defining 
the spatial geometry of the quantum computer as shown 
in Fig. [1] The number of qubits that can be directly 
connected depends on the dimensionality of the confining 
lattice. A three dimensional lattice offers the benefit of 
a higher packing density but it is also substantially more 
difficult to manipulate and measure the qubit state at a 
single site without disturbing neighboring sites in a three 
dimensional geometry We will therefore concentrate on 
a two-dimensional array of sites that can be individually 
addressed. Quite arbitrary arrays of trapping sites can 
be readily produced with diffractive optical elements or 
spatial light modulators using known techniques. We ini- 
tially assume that each site contains a single atomic qubit 
and then extend this to allow for a three dimensional sub- 
lattice within each site which is collectively addressed and 
contains an ensemble of K atoms that encode an N qubit 
register. 



Maximum number of connected qubits in a 
lattice 



A two-qubit gate between sites 1 and 2 separated by a 
distance R is achieved by exciting atoms in the two sites 
to a high lying Rydberg level with principal quantum 
number n. When the sites are separated by a relatively 
large distance R we have a van der Waals interaction 
of the form f/ v dw = Cq/R 6 . The asymptotic scaling at 
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Here q is the electronic charge, ao is the Bohr radius, 
and En is the Rydberg energy. The n scaling seen in 
Eq. ([1]) should be valid for all atomic species. However, 
this estimate turns out to be too conservative since it 
is based on assuming the Forster energy defect scales as 
6 ~ En/n 3 . In the heavy alkalis the difference between 
the s and p quantum defects is close to 0.5 which results 
in near cancelation of the 1/n 3 dependence. We find that 
numerical calculations of S and G^IIt} are reproduced 
much better for 50 < n < 200 by using S = 5E^/n 4 



which implies Cq 



This scaling will be obtained 



for Rydberg s states whenever the difference of the s and 
p quantum defects is an integer, half integer, or quar- 
ter integer. In addition, even when the quantum defects 
do not conspire to minimize the Forster energy, external 
fields can be used to Stark and/or Zeeman tune the levels 
into resonance. In light of these considerations we will 
assume that 5 = kgEn/n , with k$ a scaling constant, 
which gives 
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Although we will be interested in what follows in Ho 
atoms, due to lack of detailed knowledge of the Ho Ry- 
dberg spectrum and excited state lifetimes, we base the 
numerical estimates given in this section on Rb. 

To find the maximum number of qubits that can be 
directly connected using Rydberg interactions we first es- 
timate the maximum separation i? max for a desired gate 
error E. In the van der Waals limit we assume the exci- 
tation Rabi frequency Q is large compared to the inter- 
action frequency A v dw = U v dw/^- I n this limit, neglect- 
ing small corrections due to the finite energy separation 
of the hyperfine ground states, the minimum gate er- 
ror averaged over all two-qubit input states is[lq] E = 

q 2/3 -I 

2 1 / 3 (Ot) 2 / 3 ' w ith T the Rydberg lifetime. At room tem- 
perature the blackbody background limits the Rydberg 
lifetime to r ~ Ton 2 so the gate error scales as E ~ n -4 / 3 . 
This error is achieved when the interaction strength is 

set to the optimum value A Q p t = (f ) 1/3 = -grk^E- 
Using the van der Waals scaling the maximum qubit sep- 
aration at a fixed error is thus 
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Given i? ma x the number of sites that can be connected 
scales as (Rm&x/ D m i n ) d in d dimensions with D m - m the 
minimum usable lattice spacing. The minimum spacing 
is determined by several limits. The sites must be opti- 
cally resolvable in order to perform gate operations on a 
desired qubit without disturbing neighboring qubits. In a 
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lattice geometry there is an additional limitation related 
to the fact that we must avoid interactions between the 
highly excited electron of a Rydberg atom and a ground 
state atom at a neighboring site. This requirement can be 
written as D mm = fci<2o?i max with a safety factor k\ > 1. 
For alkali atom s states the radial wavefunction scales as 



Rn-,0 ~ e~ r / aon *U{l-n*,2,2r/a Q n*). 



with n* the effective principal quantum number and 
U the confluent hypergeometric function. At large r, 
Rn*,o ~ e -r /° ™ (r/aon*) n _1 which has a maximum at 
r n * = ao(n*) 2 . At a larger distance r' = k\r n * the wave- 
function is less than its maximum value by a factor of 
e^^-V/kf' 1 . Setting this factor to 10 2 (so the prob- 
ability density is reduced by a factor 10 4 ) and n* = 100 
we find k\ = 1.32 or D m \ n = 0.7 /im. Allowing for an 
additional positional uncertainty in the plane of the lat- 
tice for each atom of 0.15 jum suggests a lower limit of 
D min ~ 1.0 /urn, which corresponds to fci c ff = 1-89. This 
value of D min is consistent with the single site address- 
ability requirement using visible and near infrared lasers 
for internal state manipulation, provided fast diffraction 
limited addressing optics and "top-hat" shaped beams 
are used. 

With these considerations in mind we find the maxi- 
mum number of interconnected sites when relying on a 
1 /R 6 van der Waals interaction in square or cubic lattices 
is 
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We see that the number of connected qubits at a fixed 
value of the gate error scales ~ n 2 / 3 in 2D and ~ n in 3D. 
Putting in k% = 1.89, kg — 5, tq = 54 ns, a target error 
ofE = 0.001, and n = 100 gives N^ 3 ^ dW = 470, 7600 
with an array that is about 22 /mi on a side. These num- 
bers can be increased further by increasing n, but there 
are some practical limitations to letting n be arbitrarily 
large. These include the difficulty of rapid laser excita- 
tion of high lying Rydberg states as well as sensitivity 
of the Rydberg states to external fields. Excitation of 
specific levels with n as large as 500 has been achieved 
with careful shielding (26|. We will assume a conserva- 
tive limit of n m ax = 100 which ensures sufficiently long 
Rydberg radiative lifetimes that gate errors O(10 -3 ) are 
feasible (la|. 

We conclude this section by emphasizing that balanc- 
ing the requirement of strong long range interactions with 
the necessity of minimizing neighboring site ground to 



Rydberg interactions fundamentally limits the number of 
qubits that can be directly connected without mechan- 
ical motion, swap gates, or conversion to flying qubits. 
Using a set of realistic parameters and a gate error tar- 
get of E = 0.001 the limit is about 500 qubits in a 
2D geometry and 7500 qubits in a 3D geometry. While 
this number is not extremely large it may be sufficient 
for performing quantum simulations that are intractable 
on classical computers, and it significantly exceeds the 
limit of present trapped ion based approaches for which 



B. Maximum number of ensembles 

One of the primary challenges associated with building 
a 500 qubit device is the need to prepare 500 sites with 
a single atom in each. In the following we will discuss an 
approach to increasing iV max based on collective encoding 
of an N — 60 qubit register at each of a smaller number 
of sites. Collective encoding removes the requirements 
of individual qubit addressing and preparation of sites 
with single occupancy. The scaling laws are therefore 
different than those of Eqs. ((H). Let us assume that K 
atoms are used for register encoding in each ensemble 
and that the atoms are contained in a spherical volume 
of diameter Dk- We envision an architecture where gates 
can be performed within one ensemble and also between 
two ensembles in the array shown in Fig. [T] For inter- 
ensemble gates we will access Rydberg states with n — 
100 and therefore the atom spacing inside each ensemble 
must respect the limit of -D m i n = 0.7 fim found above. 

As will be discussed in Sec. |TV] we propose to base 
each ensemble on K ~ 100 atoms trapped in a lat- 
tice with periodicity Z? m i n and a filling factor / = 0.5. 
These numbers imply a sphere with diameter Dk = 

(f) 1/3 Ani„(^//) 1/3 = 5.1 fan. The relatively small 
5.1 /im maximum separation between atoms inside one 
ensemble allows intra-ensemble gates to be based on the 
blockade mode of operation which is insensitive to the 
precise value of the atom separation, while inter-ensemble 
gates can be based on the limit of A v dw <C fi discussed 
in the preceding section. An array of ensembles, each 
with Dk = 5.1 /im placed on a 2D grid with periodicity 
D = 5.3 fim gives k\ ie g = 10. and Eq. (fjh) then pre- 
dicts that 17 such sites can be connected. With each site 
containing N = 60 qubits this forms the basis for a 1000 
qubit scale processor. We will discuss in detail in Sec. 
IIVI an efficient method by which loading of the ensembles 
can be achieved in parallel. 



III. SCALING BY COLLECTIVE ENCODING 

The limit on iV max found in the preceding section can 
be increased by a factor of N using collective encoding 
of a TV qubit register in an ensemble of K > N atoms at 
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FIG. 2: (color online) Qubit encoding in the symmetric states 
of an ensemble of (2N + l)-level systems. \s) denotes the 
reservoir state. The figure depicts the state 1 10 . . . 01). 
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each lattice site. We have recently described two different 
approaches to encoding a register in an ensemble of atoms 
each with N' > N stable ground states [IJ, [HI]. In the 
approach shown in Fig. [2] the |0j) and states of each 
register qubit i are associated with single collective exci- 
tations of different hyperfine ground states. While this 
approach requires two internal states per qubit, which is 
not the most efficient encoding possible [14}, it has the 
distinct advantage that error correction protocols can be 
implemented in a fairly straightforward wayfl5j|. 

The register size is limited by the number of stable in- 
ternal states. In [l4[ we discussed the use of Cs which 
has 16 hyperfine ground states, and would allow for a 7 
qubit register using the encoding of Fig. Other atomic 
species have many more stable ground states. The rare 
earth atoms in particular with unfilled / shells have large 
nuclear and electronic spins, and therefore many hyper- 
fine ground states. The rare earths also have large hyper- 
fine splittings and ground state multiplet splittings which 
turns out to be useful for qubit preparation and readout. 
In this section we discuss the prospects for encoding a 
60 qubit register in Ho which has 128 hyperfine ground 
states, more than any other stable atom. 



A. Register encoding in hyperfine ground states of 

Ho 



FIG. 3: (color online) The known Q 111 levels of Ho below 
25,000 cm - are shown with odd parity levels indicated by 
a dashed line. Transitions suitable for cooling and trapping 
(a-g), shelving to the J = 11/2 metastable ground multiplet 
(si, S2), and readout by resonance fluorescence (n, 7-2, 7-3), are 
indicated. 

characterized by large hyperfine splittings which are con- 
venient for qubit encoding. The hyperfine constants 
of the ground state multiplets are known with high 
accuracyfl^, [3(| leading to the splittings shown in Fig. 
31 A two-state per qubit encoding giving 60 hyperfine 
qubits is shown in Fig. 2] Each bit is associated with 
the levels \F,m)(= |0»>) and \F + l,m) (= \U)) with 
F = 4, 6, 8, or 10. We see that all states, except for the 8 
with F odd and ttif — ±F are assigned to qubits. The 
energy splittings between the and 1 states of each qubit 
range from 4.3 to 8.3 GHz. The qubit phase is sensitive 
to magnetic fields since the gF Lande factors are different 
for all F levels. The values indicated in the figure assume 
L — S coupling and no configuration mixing, which is an 
accurate description for the ground state multiplet [30T] . 
A stable magnetic environment is therefore necessary to 
prevent qubit dephasing. One or more of the unused 
states can be assigned to the reservoir \s). A convenient 
choice is to use \s) = 1 1 1 , 11) this state can be populated 
by pumping on transition f) at 545.3 nm with er + light. 



To start let us recall some basic facts about the struc- 
ture of the rare earth Ho shown in Fig. [3l There is one 
stable isotope 165 Ho which has nuclear spin 1 = 7/2 and 
a ground electronic configuration 4/ 11 6s 2 . The ground 
state term is 4 I°, J = 15/2, giving 128 hyperfine states 
with 4 < F < 11. Transitions suitable for laser cool- 
ing and trapping are labeled a)-g). Ho is character- 
ized by very large term dependent shifts and as seen in 
Fig. [3] there are three odd-parity metastable levels in the 
ground multiplet with J = 13/2, 11/2, 9/2 at energies of 
5420,8605, and 10700 cm" 1 . These auxiliary levels will 
be used for qubit readout. 

The ground and low-lying excited states of Ho are 



B. Cooling and trapping 

The rare earth atoms including Ho have been magnet- 
ically trapped and collisionally cooled with He buffer gas 
to mK temperatures 31]. This approach requires strong 
magnetic fields of several T and is not directly suitable 
for qubit storage. There has also been recent progress 
in laser cooling and magneto-optical trapping of the rare 
earth Er without the need for a repumping laser, despite 
the lack of a closed cycling transition[32|]. It has been 
argued convincingly that the large magnetic moments of 
Er contribute to the success of laser cooling and trap- 
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FIG. 4: (color online) Hyperfine structure of the Ho 4/ 11 6s 2 ( 4 7i 5 /2) ground state. Qubit assignments 1-60 are indicated 
together with hyperfine splittings and g factors. 



ping, despite the fact that a cycling transition was not 
used, due to the possibility of magnetic trapping in a 
moderately strong quadrupole field. The spectrum of Ho 
is very similar to Er as concerns laser cooling possibili- 
ties, and if anything looks even more favorable due to the 
presence of sufficiently strong closed cycling transitions 
for cooling. 

The transitions labeled a)-f) in Fig. [3] are all closed 
cycling transitions between the J = 15/2 ground state 
and J = 17/2 excited states with vacuum wavelengths of 
1193,867.3,660.9,608.3,598.5, and 545.3 nm. The only 
dipole allowed decay path from the upper level of these 
transitions is back to the ground state so we expect leak- 
age out of the cycling transition to be negligible. Transi- 
tions a)-d) have an upper level configuration of 4/ 10 5cffis 2 
which is dipole allowed, but has small oscillator strengths 
to the ground state. Transitions e) and f) have an 
upper level configuration of 4/ 10 6s6p and are coupled 
more strongly with the ground state. The transition 
linewidths are however only known for d) and e) which 
havefH, [H, [H (7d,7e) = (0.25,0.92) x 10 6 s" 1 which 
give Doppler cooling limits of 0.95,3.5 /iK. These are 
attractively low temperatures but the lines may be too 
narrow to allow efficient capture from a background ther- 
mal vapor. Transition f) may have a larger linewidth, but 
this will have to be studied experimentally If not there is 
the possibility of cooling on the very strong transition la- 
beled g) to the level 4/ 11 ( 4 / 1 ° 5/2 )6s6p( 1 P 1 °), J = 17/2 at 
410.5 nm which appears analogous to the strong blue line 
used for cooling of Er to subDoppler temperatures [36| . 
This transition has a linewidth of 7 = 204 x 10 6 s" 1 
and a Doppler cooling limit of 780 /iK. The subDoppler 
temperatures observed in Er may have been due to a for- 
tuitous coincidence of ground and excited state (/-factors, 
and may not occur in Ho. Nevertheless the prospects for 
laser cooling of Ho look very promising either in a single 
step using lines d) or e), or by using g) first to have a 



large capture range, followed by deep cooling on d). 

We note that the Er cooling experiments were success- 
ful without the use of a repumper, even for 167 Er which 
has ground state hyperfine structure . The same may be 
true for Ho, although for quantum information applica- 
tions it will be necessary to prepare the sample with all 
population in the reservoir state \s). This can be done 
either by shelving the |s) level in one of the metastable 
ground multiplets, and then blowing away the unwanted 
atoms with unbalanced resonant light, or by adding re- 
pumper frequencies to deplete the lower hyperfine levels 
during cooling. With either approach it may be nec- 
essary to have as many as 7 repumping frequencies to 
depopulate all hyperfine levels 4 < F < 10. 

After laser cooling and magneto-optical trapping we 
can transfer the atoms into tightly confining optical traps 
for qubit manipulation. In order to hold ground and Ry- 
dberg state atoms in the same trap it is preferable to use 
blue detuned trapping light such that both the ground 
state and the Rydberg electron have negative polariz- 
abilities and can be simultaneously confined. A blue trap 
also minimizes spontaneous scattering of trapping light 
and is optimal for achieving long ground state coherence 
times fl8jj . Dipole allowed excitation from the Ho ground 
state to an even parity level is possible at a large number 
of wavelengths. The use of a short wavelength transition 
allows for tighter spatial focusing and confinement. 

Figure[5]shows all levels between 24000 and 2600 cm -1 , 
including the upper level of the strong g) transition dis- 
cussed above at 24360.81 cm -1 . Although this transition 
is a good candidate for cooling it turns out not to be 
directly convenient for a blue detuned optical trap due 
to the proximity of a second, almost equally strong tran- 
sition to the 4/ 11 ( 4 / 1 5/2 )6s6p( 1 P 1 °), J = 15/2 level at 

24060.80 cm -1 . This line has very similar strength to g) 
with[35[ 7 = 200 x 10 6 s _1 . In order to accurately esti- 
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FIG. 5: (color online) Levels between 24000 and 2600 cm' 1 
that are connected by dipole allowed transitions with the 
ground state are shown in solid blue, and dipole forbidden 
levels are dashed. The dipole allowed levels are labeled with 
their energy and lifetime. Part b) shows trap depth (solid 
curve) and scattering rate (dashed curve) for light of the in- 
dicated energy. 



mate trapping conditions we have accounted for all the 
levels in Fig. [5] which have dipole allowed transitions 
to the ground state. We can estimate the possible trap 
depth and photon scattering rate by assigning the indi- 
cated lifetimes to a single transition to the ground state 
for each level. This is not a bad approximation since 
there are no odd parity states between the ground mul- 
tiplets and 18572 cm -1 , so the indicated lifetimes are 
roughly equal to the inverse of the decay rate 7 to the 
electronic ground state for each level. We also neglect 
the level with J = 13/2 at 25931 cm" 1 . The lifetime is 
unknown for this level, presumably because it has con- 
figuration 4/ 10 5<i6s 2 and is only weakly coupled to the 
ground state. There are also additional dipole allowed 
transitions to levels at higher energy, the first one being 
4/ n 6s6p, J = 15/2 at 26958 cm" 1 . Again the lifetime is 
unknown, and presumably there is only weak coupling to 
the ground state. 

Using the two-level formulas for the optical potential 
U = Tf-^^-f, (w a is the transition frequency, A = lo — iv a 
is the detuning, and / is the intensity) and the scatter- 

3 2 -v 2 

ing rate r = > an d summing the contributions 



from each of the allowed levels in Fig [5^,) we obtain the 
optical potential and scattering rate curves shown in Fig. 
[5b) for 5 mW of power focused to a Gaussian waist (1/e 2 
intensity radius) of w — 5 /im. We see that the ratio 
of well depth to scattering rate is optimized at energies 
above 25500 cm" 1 . Tuning to 25700 cm" 1 (389 nm) gives 
a trap depth of 100 /iK and a spontaneous scattering rate 
of just under 10 s _1 . Note that the scattering rate will be 
further reduced in a blue detuned trap by approximately 
ksTa/U where kg is the Boltzmann constant and T a is 
the atomic temperature. We conclude that optical trap- 
ping of cold Ho looks relatively straightforward with low 
power light in the near uv part of the spectrum. 



Initialization, single qubit operations, and 
measurements 



Assuming that the task of cooling and trapping Ho has 
been completed let us now consider how a qubit register 
can be initialized and manipulated. Referring to Fig. [3] 
we see that any of the transitions a)-f) is a closed cy- 
cling transition. Using a + polarized light and including 
repumping frequencies to depopulate all hyperfinc levels 
with F < 11 will result in initialization of all atoms into 
the reservoir state \s) = 1 11, 11) . The Ho ground state has 
vanishing configurational mixing with other levels [13] so 
the limiting factor as regards the efficacy of state prepa- 
ration will be due to mixing of the upper state of the 
transition used for pumping, as well as Raman transi- 
tions into lower hyperfine states. As mentioned above, 
one approach to maximizing the state purity is to shelve 
the \s) state in an excited metastable level of the ground 
electronic configuration using transitions s\, S2 in Fig. [3] 
at 586.2 and 1183 nm respectively, which proceed via the 
4/ n ( 4 J° 5/2 )6s6p( 3 P 1 °), J = 13/2 excited state. We can 
then push away any unwanted population left behind in 
the other ground state levels using unbalanced force from 
light resonant with the e) transition. 

Having prepared all atoms in |s) the qubit register 
must be initialized to a fiducial state. This can be done 
using sequences of stimulated Raman transitions start- 
ing with the state "furthest" away from |11, 11) and then 
working backwards. Since the register states are collec- 
tive states the initial step must involve a Rydberg inter- 
action, after which additional steps can be taken using 
stimulated Raman via for example transition e). Specif- 
ically, starting in | s) we prepare the collective state with 
unit excitation in |Oeo) = 1 10, 10) via Rydberg blockade. 
This state can then be mapped onto |0i) = |4, —4) using 
7 steps of two-photon stimulated Raman via the cycling 
transition e). We then prepare | O2) , IO3), ... and so on. 

Single qubit rotations |0i) <-> |li) on register bit i 
are then straightforwardly performed using two-photon 
stimulated Raman beams tuned close to transition e). 
The qubit state energy separations range from 4.3 GHz 
for bits 1-9 (F = 4 -> 5), to 8.3 GHz for bits 40-60 
(F = 10— > 11). In order to distinguish between different 
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bits with the same value of F we rely on Zeeman shifts 
due to an external magnetic field. Since the gp factors 
are different for each F level this is possible. 

Qubit readout can be performed using shelving. We 
transfer the part of bit i to the metastable J = 11/2 
ground state multiplet using two-photon stimulated Ra- 
man with fields si, sq. as described above, the only change 
being that the frequency "si — S2" must be tuned to be 
resonant with the transition from level to a level in 
the J = 11/2 multiplet. The J = 11/2 level can then 
be coupled to 4/ 11 ( 4 J° 5 y 2 )6s6p( 3 P2 ) using the n transi- 
tion at 878.8 nm to generate fluorescence photons. The 
upper level only has dipole allowed decay paths to the 
metastable J = 9/2,11/2,13/2 levels shown in Fig. [3j 
all of which can be pumped back up to the same up- 
per level. This requires two additional repumpers r2,r% 
at 1231 and 746.2 nm. If the measurement gives a null 
result the register bit is projected into |0j) and we are 
done. 

If the measurement result is |lj) then the register bit 
will have to be restored for further processing. The atom 
can be optically pumped into the | J = 11/2, F = 9, mj? = 
9) state using <r + light on transition r\ and then coher- 
ently transferred back to \s) using transitions s\ and S2- 
All bits j > i can then be swapped down to fill the hole 
at bit i flOi+i) \0i), -> |T), |0 i+2 ) -> |0 i+ i), 

lli+2) — ► etc.), and the mapping between infor- 

mation and register bits can be relabeled to account for 
the change. Bit 60 can then be reinitialized from |s) via 
a blockade operation and the computation can continue. 
Alternatively, all the swap operations can be avoided by 
restoring the atom to the reservoir \s) and then per- 
forming the first blockade mediated operation from \s) 
back to a state in J = 11/2. From there the state can 
be coherently transferred to \F\ mjf) in the J = 11/2 
level using stimulated Raman. The measured register 
bit I J = 15/2, F, mp) can thus be restored after the mea- 
surement to I J = 11/2, F',m F >) where \F' -F\<2 and 
\mpi — mp\ < 2. This state can then be returned to re- 
store the register bit | J = 15/2, F,mp) using stimulated 
Raman on si,S2- This eliminates the need for swapping 
many register bits at the expense of an additional laser 
to couple J = 11/2 to Rydberg levels. Using either ap- 
proach we thus have a protocol which allows for resetting 
of a measured qubit so that a computation can proceed 
as long as there is no physical loss of atoms. 

Finally, there are two additional issues which should 
be mentioned in connection with single qubit operations 
and measurements. The first is that the optical trapping 
discussed above using light at 389 nm will not produce 
the same trapping potential for the J = 9/2, 11/2, 13/2 
states. This can potentially be solved by adding addi- 
tional beams that are each tuned to create an attractive 
potential for these states, but are far off-resonant with 
respect to the ground state, and do not disturb the un- 
measured part of the register. Since these beams are only 
needed during a brief measurement some amount of pho- 
ton scattering can be tolerated in J = 9/2,11/2,13/2, 



and the choice of wavelengths will not be as constrained 
as for the calculations shown in Fig. O It is also feasible 
to consider readout based on cross entanglement between 
atoms of different species [13], which removes the require- 
ment of using the metastable levels entirely. A detailed 
study of this possibility is outside the scope of this work. 

The second issue is related to Zeeman selectivity. The 
ground state hyperfine separations between levels with 
neighboring values of F range from 4.31 to 8.28 GHz, so 
selective shelving to or from the auxiliary levels as well 
as Rydberg levels are well resolved as regards different 
F values. For each value of F, bits with different m 
can be isolated by applying a magnetic field B to give 
shifts AUip jm \ — gpfisBm. The gp values range from 
0.82 for F = 11 to 2.04 for F = 4. Moderate fields of 
under 2.2 Gauss will thus be sufficient to get 2.5 MHz 
of separation between all qubits as regards shelving and 
Rydberg excitation. This allows the use of greater than 
100 kHz Rabi frequencies while keeping the probability of 
population transfer of a nonaddressed bit below the 10~ 3 
level in the worst case, and much less for the majority of 
register bits. 

There remains, however, the problem of selectively ad- 
dressing bits for single qubit rotations. In this case the 
resonance condition scales as AU\F, m ) — AC7|^_ l jrl ^ — 
(gp — gp-ijUBBm. The worst case is for bits 40-60 since 
<7io = 0.86 and <?n = 0.82. The selectivity between 
bits with neighboring m values is thus only 0.044/isB 
or about 61 kHz/G. To run single qubit operations at 
a 100 kHz Rabi frequency and have errors at the 10~ 3 
level on a neighboring bit requires a detuning of about 
3 MHz or about a 50 G field. While this is not particu- 
larly large we simultaneously require small dephasing on 
the most sensitive bits 1 and 9 which see a big differen- 
tial shift of 135 MHz. Any known differential rotation can 
be accounted for, so the feasibility of Zeeman selectivity 
without inducing unwanted dephasing rests on the abil- 
ity to create a very low noise bias magnetic field. Such a 
large field only needs to be turned on for a time given by 
the inverse Rabi frequency or about 10 (is. To keep the 
dephasing error at 10 -3 of a radian would require a field 
stability of - 10~ 3 /1350 - 10~ 6 . This is not impossible 
but will be a technically challenging requirement. 

There are several possible approaches to mitigating 
this stability requirement. Since shelving is well resolved 
by much smaller magnetic fields than those needed to re- 
solve single qubit rotations, we can shelve bits \0i), |T) 
to the metastable J = 11/2 level, do the rotation there, 
and then return them to the ground state levels. Al- 
ternatively, to perform a rotation on a qubit encoded in 
F = 10, 11 states we first swap the bit with a bit en- 
coded in F = 8,9 (11 -> 9 and 10 -> 8), relabel the 
information, and then perform the rotation on the bit in 
F = 8,9. This has the advantage that the swap opera- 
tion has a larger selectivity governed by the gp of the 
lower F states. Swapping in this way from F = 10, 11 to 
F = 8, 9 reduces the bias field requirement by a factor of 
(59-5n)/(5io -511) = 2.3. 
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D. Rydberg interactions 

Finally we need to consider Rydberg interactions in 
Ho. The asymptotic scaling of Eq. (J2]) is expected to 
hold at large n for any atomic species with singly excited 
Rydberg states so the estimates found in Sec. [TT1 should 
remain valid, although the precise numerical values will 
require adjustment. We know of no principal reasons why 
Rydberg blockade analogous to what has been studied in 
the alkalis [23j should not be possible but a detailed char- 
acterization remains a topic for future investigation. We 
are aware of only one experimental study of the Rydberg 
structure of Ho[38j where the Rydberg series correspond- 
ing to excitation of 4/ 11 6snp was observed by collisional 
ionization. Resolved Rydberg levels up to n ~ 47 were 
seen. As these are even parity states excitation from 
the ground state requires either one ~ 210 nm photon 
or a three step excitation at longer wavelengths. The 
4f 11 6sns or 4/ 11 6snd series which are accessible by two 
step excitation are a second possibility. 

The largest uncertainty as regards the feasibility of Ry- 
dberg gates in Ho concerns the Rydberg state lifetimes. 
In Sec. |TT] we assumed r ~ n 2 blackbody scaling. If the 
Ho Rydberg series are perturbed by interactions between 
the valence and core electrons the lifetime could be sub- 
stantially altered. If strong core - valence interactions do 
occur they are likely to be specific to particular n val- 
ues, so that with judicious choice of the Rydberg level 
it should be possible to minimize the impact of series 
perturbations. 



IV. PREPARATION OF A LATTICE OF 
COLLISION SUPPRESSED ENSEMBLES 

In light of the principle feasibility of collective encoding 
in Ho discussed above, it is interesting to examine how 
an array of ensembles can be efficiently prepared in the 
geometry shown in Fig. [T] 

The use of a many atom ensemble at each site, in- 
stead of requiring deterministic loading of single atoms, 
immediately removes one of the prime challenges of neu- 
tral atom quantum computing, which has been the dif- 
ficulty of preparing a singly occupied lattice of optically 
resolvable sites. There has been a considerable amount 
of work on this problem based on several different solu- 
tions. An elegant approach relies on the BEC to Mott 
insulator transiti on p3| as has been demonstrated in re- 
cent experiments [40j. Unfortunately the transition only 
works in finite time with a very short period lattice 
which does not give optically resolvable lattice sites. This 
has been partially addressed by the demonstrated trans- 
fer of the insulating state to a longer period lattice [4l|. 
While the Mott insulator transition at finite temperature 
does not result in a lattice that is filled with perfect fi- 
delity there are also possible solutions to purifying the 
lattice[!2]. These purifying steps can also potentially be 
used on a lattice which is imperfectly filled directly from 



a ma gnet o-optical trap (MOT), without needing a BEC 
phase |43j. One of us has also proposed a determinis- 
tic loading scheme[44{ which relies on Rydberg blockade 
to remove all but one atom from a multiply filled site. 
This can potentially be implemented in parallel on a large 
number of sites. 

The ensemble based approach described here is poten- 
tially much simpler since it does not require preparing 
sites with single atom occupancy. Based on the results 
of[14j a 60 bit register, including error correction, needs 
an ensemble of K > 60 atoms. For the purposes of the 
present discussion we will target a value somewhat larger 
say, K ~ 100. These K atoms should all be confined to a 
small volume to allow an effective Rydberg blockade and 
to maximize the total number of processor qubits. The 
disadvantage of a multi-atom ensemble is that high den- 
sity samples suitable for Rydberg blockade experiments 
are susceptible to collisions which will drastically reduce 
coherence times. To solve this problem we propose to 
load the ensemble in a short period lattice, thereby elim- 
inating collisional decoherence. 

The loading protocol is shown in Fig. O We start with 
Ho in a MOT, and assume Doppler cooling to a few f/K 
using transition e). The atoms are then transferred into 
a 3D blue detuned lattice (3 pairs of beams with orthogo- 
nal polarizations, and frequency shifts to avoid unwanted 
interference) created with A = 389 nm light, as described 
above. The angle between each pair of beams is adjusted 
to give a lattice spacing of A = 0.7 fim which corresponds 
to a lattice site density of notice = 2.9 x 10 12 cm" 3 . 
Magneto-optical trapping of Er has resulted in peak den- 
sities of n a ~ 2 x 10 9 cm -3 and we expect similar base- 
line results for Ho. It is unknown what densities can 
be achieved using transient magnetic compression and/or 
evaporative cooling, but by analogy with experience from 
the alkalis where evaporation has produced small samples 
with densities above 10 15 cm _3 [45| preparing a few thou- 
sand atoms with n a ~ 1.5 x 10 12 cm -3 appears realistic. 
Loading such a sample into the above lattice would give 
a filling fraction of n a /notice ~ 0.5. Assuming Poisso- 
nian loading statistics the probability of double site oc- 
cupancy at this filling fraction is 0.09 which implies that 
for K ~ 100 about 18 sites will have more than one atom. 
The atoms in these sites will rapidly redistribute or be 
lost from the trap due to hyperfine changing collisions 
leaving a sample with slightly smaller K but only single 
atom occupancy at all sites. 

We then superimpose a blue-detuned bottle beam 
trap [46[ and drop the lattice to remove atoms outside the 
region defined by the bottle beam (frames b,c in Fig. 
The lattice is then restored and the bottle beam turned 
off (frames d,e) which leaves us with K atoms in the lat- 
tice, in a well defined spatial region. We can achieve an 
average value of K = 100 using n a ~ 1.5 x 10 12 cm 3 
and a spherical bottle beam trapping region with diam- 
eter D = 5.1 /im (the bottle beam trapping potential 
is not far from spherical for such small trap volumes). 
The probability of doubly occupied sites, and hence col- 
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FIG. 6: Protocol for loading small ensembles into a patterned lattice region. The small open circles are 
repulsive lattice sites, and the large green circle is the outline of the bottle-beam trap. The lattice is 
drawn to scale for clarity. See text for details. 
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lisions, will be somewhat higher than in the initial lattice 
state since in the bottle beam transfer phase when the 
lattice is dropped the atoms will tend to fall towards the 
center and compress. This effect will be smaller than in 
a harmonic trap since our bottle beam design provides 



a quartic transverse potential, U 



which provides 



only weak radial compression. The potential is harmonic 
axially, but by judicious choice of the time window for 
the bottle beam transfer At ~ pi V ib, where i V ib is the 
axial vibrational period, and p is an integer multiplier, 
compression effects should be minimized. 

Following the above procedure we can prepare a lo- 
calized region with K ~ 100 atoms. This can readily 
be done in parallel at multiple sites by superimposing 
an array of bottle beams on the 3D lattice. Poissonian 
loading statistics will imply ~ 10% variations in K from 
site to site, and hence variations in Rabi frequencies for 
qubit operations involving Rydberg levels. For the first 
encoded qubit which is supported by K ~ 100 atoms 
the Rabi frequency variations will be \J\.\ ~ 1.05 or 
±5%. For the last encoded qubit which is supported 
by K — 59 ~ 40 atoms the Rabi frequency variations 
will be vT25 ~ 1.12 or ±12%. These variation will not 
prevent high-accuracy gate operations since composite 
pulse sequences can be readily used to remove the de- 
pendency on Rabi frequency variations. We have thus 
arrived at the situation shown in Fig. [T] with an array 
of closely spaced ensembles. If the diameter of each en- 
semble is about 5.1 /urn then we might space them in a 
plane by D ~ 5.3 /im using addressing beams with flatter 
than Gaussian profiles to minimize site to site cross talk. 
The scaling arguments of Sec. [TT] imply we can connect 
roughly 17 ensembles. The total processor size is thus 
60 x 17 = 1020 directly coupled qubits. 



V. DISCUSSION 

We have presented a path towards a neutral atom gate 
array that could allow as many as a thousand qubits to be 
interconnected without resorting to mechanical motion, 
swap chains, or interconversion between stationary and 
flying qubits. The scaling relies on optimization of qubit 



array and Rydberg parameters to allow up to 500 sites to 
be directly coupled in a planar geometry. The estimate 
was based on realistic parameters, including gate errors 
of E = 0.001 and a maximum principal quantum num- 
ber of n = 100. This limit may be too conservative, and 
the ultimate limit in a planar geometry may be higher 
than predicted here. Nevertheless manipulating single 
atoms in such a large number of sites presents a difficult 
experimental challenge. We can avoid the overhead asso- 
ciated with individual addressing and single atom load- 
ing at each site by invoking collective encoding. Using 
Ho atoms this gives a multiplying factor of 60 per site. 
Since the ensemble sites are larger than those holding a 
single atom we arrive at a limit of about 17 sites, each 
with 60 collectively addressed qubits, giving a 1000 qubit 
scale device. Our estimates are based on available spec- 
troscopic data for Ho. More information, particularly 
concerning the Rydberg states of Ho, will be needed to 
further refine the limits discussed here. 

Another attractive application of collective encoding 
would be to use a single ensemble with a 60 qubit regis- 
ter as a quantum repeater, or as part of a hybrid quantum 
computation scheme involving small, optically connected 
registers fioj. The 60 qubits are sufficient for 8 logical 
qubits, each built from 7 physical qubits, or, for exam- 
ple, 5 logical qubits, each built from 5 decoherence free 
subspace encoded qubits, that in turn require 2 physi- 
cal qubits each 47) • In this way 50 physical qubits would 
enable 5 error corrected and potentially low decoherence 
logical qubits in one ensemble. The ensemble qubits can 
be efficiently mapped to and from photonic bits using 
blockade mediated preparation of excited states [44(, and 
the large choice of transitions in Ho provides flexible op- 
portunities for coupling to wavelengths compatible with 
optical fiber transmission. 

There are undeniably many challenges associated with 
collective encoding, particularly in rare earth atoms that 
are relatively poorly studied and have not been widely 
used for laser cooling. A large number of lasers of dif- 
ferent wavelengths and frequencies are required for the 
various internal state manipulations. In some sense we 
have transferred the complexity of moving the informa- 
tion about spatially, as in current approaches to scalable 
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ion traps, to the problem of dexterously moving informa- 
tion between internal states. Although the overall com- 
plexity required to build a scalable quantum processor 
will remain very high, we believe it is worthwhile to ex- 
plore a wide range of approaches. Indeed, we are not 



aware of any approaches to building a thousand qubit 
scale quantum logic device that are simple. 

This work was supported by the NSF and ARO- 
IARPA. M. S. thanks Jim Lawler for helpful discussions 
on Ho spectroscopy. 
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